Disorder Induced Ferromagnetism in Restricted Geometries 
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We study the influence of on-site disorder on the magnetic properties of the ground state of the 
infinite U Hubbard model for restricted geometries. We find that for two dimensional systems disor- 
der enhances the spin polarization of the system . The tendency of disorder to enhance magnetism 
in the ground state may be relevant to recent experimental observations of spin polarized ground 
states in quantum dots and small metallic grains. 
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The interplay between disorder and interactions |l]] 
and the possibility that it leads to ground state ferro- 
magnetism has been the subject of much interest In 
several new experiments in restricted geometries, such as 
zero temperature transport measurements of the conduc- 
tance through semiconducting quantum dots || and car- 
bon nanotubes 0] , tantalizing hints of a weakly ferromag- 
netic ground state of small systems with a few hundreds 
of electrons have appeared. The ground state spin po- 
larization may be directly measured by coupling the dot 
or tube to external leads and measuring the differential 
conductance. Also from a recent mean field treatment 
of electron-electron interactions in disordered electronic 
systems [5] as well as from a numerical study of such 
systems [6] a partially magnetized ground state seems 
probable. 

The canonical model for the study of itinerant fer- 
romagnetism is the Hubbard model, described by the 
Hamiltonian 
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T a ja + c.c. + U m 
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where a ia is the fermionic creation operator on site i with 
spin <t, riij = a\^a^, and the on-site energies £j are drawn 
randomly according to a uniform distribution between 
— W/2 and W/2. The large U regime of the model has 
attracted much interest due to its relevance to the the- 
ory of high-T c superconductivity . Although the model 
clearly does not contain many of the physical attributes 
of the typical experimental system such as a quantum 
dot (especially at the infinite U limit), nevertheless, it is 
important to gain insight into the complicated problem 
of the influence of disorder on the spin structure of in- 
teracting electrons in restricted geometries by studying 
simplified models. Moreover, the infinite U limit has the 
attractive feature of suppressing antiferromagnetic corre- 
lations which are clearly not relevant to quantum dots, 
even in the clean limit j^]. 

In this paper we study the effect of disorder on the en- 
hancement of the spin polarizability of small interacting 
systems. For one dimensional infinite U Hubbard mod- 



els, disorder does not essentially change the magnetiza- 
tion behavior j^] . On the other hand, for two dimensional 
systems, it is found that sufficiently strong disorder sup- 
presses the singlet-favoring effects and the spin polariza- 
tion of the system increases. From analytical arguments 
and numerical calculations it is shown that in the pres- 
ence of disorder the fully polarized regime extends to 
higher densities of holes. Even beyond the fully polar- 
ized regime disorder creates a tendency towards non-zero 
magnetic moments in the ground state, which is consis- 
tent with other indications for such a behavior in disor- 
dered interacting systems. 

The extensively studied Hubbard model, which is the 
simplest model of strongly correlated electrons, was orig- 
inally introduced to explain ferromagnetism |LT| |. How- 
ever, till now, little is known about the phase diagram 
of the model even at T = 0. There are a few rigorous 
results, mostly restricted to the one-dimensional model, 
or to the half-filled case. Lieb and Mattis have proven 
that for one dimensional systems (ID) for even number 
of electrons, with interaction strength U < oo and open 
boundary conditions (BC) (l2) the ground state (GS) is a 
singlet. For higher dimensions and half filling it has been 
shown 1 13 1 that for a bipartite lattice with Na {Nb) sites 
on sub-lattice A(B), the GS is non-degenerate and has to- 
tal spin \\N A — N B \, as long as U > 0. In particular, 
for the square lattice, the GS is a singlet. In the large-!/ 
limit the problem is mapped onto the Heisenberg Hamil- 
tonian, leading to antifcrromagnetism (AFM) with long 
range order ]i4| . 

An important milestone in the research of ferromag- 
netism in Hubbard models is the work of Nagaoka [15| . 
It showed that for most lattices, with nearest-neighbor 
hopping and infinite on-site interaction U, the GS is the 
fully saturated ferromagnetic state, for the case of one 
hole in an otherwise half-filled band. An extensive work 
was done in order to find whether this result can be ex- 
tended to higher hole density, or to finite U. It was shown 
that the two-hole case the GS is a singlet (l^Jl^l , but this 
GS is degenerate with the Nagaoka state in the thermo- 
dynamic limit. Various variational wavefunctions were 
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suggested to test the stability of the Nagaoka state (See 
JTs] and references therein). Bounds were given to the 
holes density for which stability may remain. The best 
bound to date is S cr < 0.2514 |2(J (where 5 is the number 
of holes per site). Still, the stability of the Nagaoka state, 
and the possibility of explaining ferromagnetism by it, is 
an unresolved problem p3| . 

In this letter we wish to show that in some sense the 
situation in the disordered case is simpler. Let us first 
sketch the situation in ID. If periodic BC are imposed in 
ID, the problem of m interacting electrons (at U = oo) 
can be mapped onto a system of m non-interacting spin- 
less fermions on ID ring with fluxes $/$o = 2-7T j/m 
(j = 0, , m — 1), where the GS energy corresponds to the 
flux j with the lowest energy ||. This mapping shows 
that the effect of the spin background in ID is trivial, 
and does not depend of the strength of disorder. 

One might expect the effects of disorder in 2D to be 
smaller than in ID, and thus no influence of disorder in 
2D as well. However, the insensitivity of the ID GS spin 
structure to disorder is accounted for by the fact that the 
spin permutation subgroup induced by the ID hopping 
terms is cyclic ||. Hence, spin background effects in ID 
are not major. On the other hand, the permutation sub- 
group induced by the 2D hopping terms is non-abelian, 
and therefore the spin background has non trivial effects 
on the dynamics of the holes. Thus we might expect an 
interplay between disorder and the behavior of the spin 
background in 2D. 

We start by considering the influence of disorder on the 
one hole case. The hopping of the hole around the lattice 
induces permutations in the spin ordering. The hopping 
term is then effectively reduced by a factor proportional 
to the expectation value of the different permutations. In 
order to minimize kinetic energy, this overlap should be 
maximal, and this is achieved by the fully polarized state 
for which the spin wavefunction is unchanged by permu- 
tations of different spins. This argument does not change 
due to disorder, and it leads to the Nagaoka's theorem 
which assures us that the GS is fully polarized, even in 
the presence of disorder. 

The two holes case is much more complicated. Al- 
though the above argument for preferring the FM or- 
der applies for the case of several holes equally, it is 
known that in the ordered case the GS is a singlet Jl6],|l7j . 
This is accounted for by the re-ordering of the spin back- 
ground in order to mask the fermionic BC between the 
holes. Thus, the spin background takes care for the 
anti-symmetrization of the many body (in fact, many- 
holes) wave function, and thus the spatial function has 
less nodes, which decrease its energy. It was shown, for 
a special variational wave function, that the resulting 
energy gain supersedes the energy increase at the bot- 
tom of the band, coming from the reduction of the hop- 
ping amplitude due to the Nagaoka effect jl6| . Thus, the 
tendency towards ferromagnetism is suppressed, and the 



spin structure, if any, is of a much more complex form 

m 

As disorder increases, the single particle functions be- 
come more and more localized (in the participation ratio 
sense). The overlap of the different single particle func- 
tions decreases, and thus, the fermionic BC constraint 
(i.e., a zero of the many body wavefunction wherever two 
particles are on the same site) becomes less restrictive, 
and does not change the many body energy much. There- 
fore, one may expect that the incentive for re-ordering of 
the spin background decreases, while, as in the one hole 
case, there still is a contribution from the hopping am- 
plitude leading to a Nagaoka state. 

Exact diagonalization for the full many-particle Hamil- 
tonian of Eq. (0) was used to test the above arguments. 
Although we have used small systems one may expect 
that due to the chaotic nature of the dots || the depen- 
dence on the number of electrons or the BC will play a 
less important role for disordered systems than in clean 
ones JhJ. Thus, the study of a small number of elec- 
trons is still useful in understanding the properties of 
dots which are populated by an order of magnitude more 
electrons. We have used up to 14 electrons on up to 4 x 4 
lattices. The size of the Hilbert space is then 471435600, 
which is far beyond exact diagonalization capabilities. 
Fortunately one can omit the double occupied states for 
U = oo and use the spin symmetry of the Hamiltonian 
to reduce this number considerably. The number of spa- 
tial functions in this case is 120, and the number of total 
spin configurations in the S z = sector is 3432, yield- 
ing a total of 411840 states. We have used group theory 
to construct the definite S states, and to decompose the 
space into subspaces of definite S and S z . The largest 
sectors (S = 1,2) consisted of 1001 spin functions and a 
total of 120120 basis functions. Group theory was used 
for constructing the matrices describing the effect of hop- 
ping on the different spin functions. We then employed 
the Lanczos algorithm to find the exact GS for 600 re- 
alizations at every disorder value. In the ordered case, 
the GS was a singlet, in accordance with |l(|[l7j]. Figure 
|l| presents the GS-spin distributions as a function of W, 
for 14 electrons on a hard- wall 4x4 lattice. The average 
spin (S) is also plotted against W, and one can see that it 
increases significantly with W . In the presence of disor- 
der, one gets a distribution of GS-spin values. For weak 
disorder, the main effect is smearing the peak at S — to 
low S values. Thus, a tendency towards weak ferromag- 
netism is clearly demonstrated even for weak disorder 
(W — 3t) which corresponds to a ballistic (mean free 
path larger than the system size) regime. Moreover, as 
disorder increases, high S values dominate the distribu- 
tion. For W — 6t corresponding to a diffusive regime a 
clear dominance of the high spin state appears. 

Similar behavior was obtained for smaller lattices and 
periodic BC. Fig. || presents the results for the same 
conditions as in Fig. R] employing periodic BC. Clearly, 
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the tendency towards ferromagnetic behavior persists, al- 
though higher values of W needed to obtain similar values 
of spin polarization. This is the result of the fact that 
for periodic BC, higher values of W are needed to gen- 
erate the same value of dimcnsionless conductance. One 
sees that, in contrast with the situation in the ordered 
case, our results are not sensitive to the lattice size or the 
BC. This manifests the chaotic nature of the dot, which 
suppress dependencies on the details of the system. 



is of order one, both distributions are quite similar, where 
(S) = 5.90 for hard-wall BC and (S) = 3.96 for periodic 
BC. 

Exact diagonalization also confirms the tendency to- 
wards non-zero ground state spin values even for a higher 
number of holes. In Fig. ^ we depict the spin distribution 
for 12 electrons on a hard- wall 5x3 lattice (3 holes). The 
GS-spin is significantly enhanced as function of disorder, 
although the most probable spin state is not fully ferro- 
magnetic. This tendency towards partial polarization of 
the ground state persists in higher hole ratios. 
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FIG. 1. The spin distributions as a function of disorder W 
for a 4 x 4 lattice with 14 electrons. For each W , the bar chart 
represents the probability of finding the GS of the system at 
a particular value of S. The inset presents the average spin 
(S) as a function of W. 
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FIG. 2. Same as Fig. | for periodic BC. 

A clear manifestation of this point is presented by the 
results for 13 electrons on a 5 x 3 lattice. In the ordered 
case, the behavior of this cluster depends dramatically 
on the BC. For hard-wall BC, the GS is fully polarized 
(i.e., S = 13/2), while for periodic BC, the GS has the 
minimal spin S = 1/2. On the other hand, once the 
system is diffusive the GS-spin polarization distributions 
become closer and when the dimensionless conductance 
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FIG. 3. The spin distributions as a function of disorder W 
for a 5 x 3 lattice with 12 electrons. 

The method of exact diagonalization is restricted to 
small lattices. In order to learn whether the tendency to- 
wards ferromagnetism persists for larger systems we turn 
to a variational method. Many authors have considered 
various variational wave functions to study the instability 
of the Nagaoka state of the U = oo model for a thermody- 
namic concentration of holes Jl8|-|2^1 . Since the reliability 
of these functions for an accurate calculation of the phase 
boundary of ferromagnetism is doubtful, we only use this 
method to get a hint about disorder influence of the sta- 
bility. For this purpose, we use the most simple of these 
functions , which is one of a single particle excitation. 
An up spin electron is removed from the occupied states 
and placed with flipped spin into another state. Direct 
calculation of the excitation energy in the ordered case 
pif yields stability of the Nagaoka state with respect to 
spin flip for 6 — 0.49 (for a square lattice), while for 
smaller hole concentration, the Nagaoka state remains 
stable with respect to this excitation. We have done the 
calculation for the disordered case by taking different re- 
alizations of disorder of a 24 x 24 system, diagonalizing 
the single particle (non-interacting) Hamiltonian to find 
its eigenvalues and eigenvectors, and then calculating di- 
rectly the excitation energy of the single flip variational 
wave function. Figure ^ shows the stability regime in 
the 5 — W plane, as follows from this excitation calcula- 
tion. For an ordered system we get the result of Rcf. Ill] 
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that the ferromagnetic state is stable for 6 < 0.49. How- 
ever, as disorder increases, the stability regime grows. It 
therefore seems that the exact results for small systems 
characterize the behavior in larger systems as well. We 
note that analytical calculation of the excitation energy 
in the disordered case according to random vector model 
(RVM) yields a completely different behavior. This stems 
from the fact that RVM ignores correlations between the 
wavefunctions and the eigenvalues which arc important 
in this case. 




Stable FM (?) 



FIG. 4. Stability curve for the single flip excitation: the 
critical hole density S vs. the disorder distribution width W 

We would like to add a remark about the U = oo limit. 
The Nagaoka's effect and the decrease in the singlet fa- 
voring effect described here, are not unique to the U = oo 
limit. However, for the Hubbard model, ferromagnetism 
arises (even for one hole) only for U S> t. The reason 
is that due to the perfect nesting property of the lat- 
tice model, the GS of the almost half filled case tends to 
be AFM. In order to wash out this tendency, the limit 
U = oo is taken. However, in reality, quantum dots do 
not show AFM behavior, since they are not described by 
a perfect lattice. One then might expect the previously 
described effect, namely the formation of larger magnetic 
moments due to disorder, to show in real quantum dots 
even for moderate values of U. 

In conclusion, the influence of disorder on the magnetic 
properties of the GS was studied. For an ordered sys- 
tem, large magnetic moments are generally suppressed, 
and the spin structure of the GS, if any, is very compli- 
cated. On the other hand, we have shown that disorder 
plays an important role in determining the spin polar- 
ization of 2D systems described by the infinite U Hub- 
bard model. Weak disorder tends to create a partially 
polarized ground state, while stronger disorder tends to 
stabilize a fully ferromagnetic GS. This behavior clearly 
indicates that there is a basis to expect that for more re- 
alistic descriptions of the experimental systems (U oo) 
disorder will play an important role in creating a spin 
polarized ground state. 
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